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£NJ , Abstract 

Strong correlations between the number of emitted jets and their average trans- 

^ \ verse momentum are found for the events resulting from the exchange of a single BFKL 

^^O . pomeron. 
00 

As well-known, strong correlations are observed experimentally between the average pt 
and multiplicities of particles produced in high-energy hadronic collisions [1]. Average pt 
grows with multiplicity. This fact can be interpreted as a consequence of multiple hard 
collisions, which result in a larger number of particles produced and a broadening of the pt 
spectrum [2]. An alternative interpretation can be given in terms of colour strings which 
i" 1 ■ are stretched between the colliding hadrons during the collision. Here larger multiplicities 

correspond to a larger number of strings, which again leads to a broadening of the px spectrum 
either because of the accumulation of more transverse momentum from the parent partons or 
because of the interaction between strings [3]. In both cases it is tacitly assumed that with 
only one hard collision or, alternatively, with only one pair of non-interacting strings there are 
no correlations between < pt > and multiplicity. Theoretically this assumption can only be 
tested within the BFKL dynamics, which presents a detailed description of particle (actually 
jet) production at high energies under certain simplifying assumptions (a fixed small coupling 
constant). The present calculation is aimed to see whether there exist correlations between 
< pt > and the number of produced jets in the hard pomeron described by the BFKL chain 
of interacting reggeized gluons. 

Our study is closely related to the paper by J.Kwiecinski, C. A. M.Lewis and A.D.Martin, 
who calculated the exclusive probabilities to observe a given number of jets from the ex- 
changed hard pomeron [4]. We shall extensively use both their method of calculations and 
some of their parametrizations. 

The results we present in this letter demonstrate that, in fact, at realistic energies in the 
BFKL chain there are strong correlations between the number n of produced jets and the 
average < pt > n for these events. We find that < px > n grows roughly linearly with n 
with a slope independent of Q 2 in deep inelastic scattering (DIS) and of the same value for 
purely hadronic collisions. The slope diminishes with rapidity y. So in the limit y — ► oo the 
correlations are expected to disappear. However, at y = 15 (x ~ 1CT 7 ) the value of the slope 
is still ~ 0.8 GeV/c. These results leave open the question of correlations between < p^ > n 
and n in the soft pomeron and in the colour string at the currently available energies. 

1 On leave of absence from St. Petersburg State University (Russia) 
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We start by recalling the main points of the formalism employed in [4] to calculate the 
exclusive probabilities for the production of a given number of jets. To facilitate comparison 
with [4] we shall borrow their notations. Let the amputated BFKL amplitude be f(y, k), when 
y is the rapidity and k is the 2-dimensional transverse momentum of the virtual (reggeized) 
gluon. Function f(y,k)/k 2 is interpreted as an unintegrated gluon distribution, related to 
the conventional gluon distribution by 

xG(x,Q 2 ) = j ^L/(In-,fc). (1) 
The BFKL equation for / may be written in the form 

f(y,k) = f<°\y,k)+a s j o y dy 1 J ^(^f(y lM ) - f(y, k )0{k 2 - q 2 )). (2) 

Here a s = 3a s /TT, q = k — k\ is the transverse momentum of the emitted (real) gluon and 
it was assumed in [4] that the driving term /(°) (the impact factor of the target) may also 
depend on rapidity. To suppress the physically unknown infrared domain and make the 
equation numerically tractable, the integration over k\ was constrained in [4] to the interval 

Ql < kf < Qj, (3) 

with Qo = 1 GeV/c and Qf = 100 GeV/c. We shall also impose this constraint, implicit in 
the following equations. 

Defining as an observable jet a real gluon with q 2 > fi 2 , one splits the integration over 
momenta in (2) into two parts by introducing 

e{q 2 -^ 2 )+e{ l i 2 -q 2 ) = l (4) 

inside the integral. The whole integration kernel is thus split into two parts: a resolved one, 
Kr, corresponding to emitted gluons with q 2 > fi 2 , and an unresolved one, Kuv-, which 
combines emission of gluons with q 2 < fi 2 and the subtraction term in (2). Explicitly, the 
action in the momentum space of the two kernels is described by 

(K R f) (k) = a s k 2 J ^^0(q 2 - fi 2 )f(ki), (5) 

(K UV f) (k) = a s k 2 J (V - q 2 )f( kl ) - f 2 6{k 2 - q 2 )f(k)) . (6) 

Exclusive probabilities to produce n jets are obviously obtained by introducing n operators 
Kr between the Green functions of the BFKL equations with kernel K\jy. If one presents 
the full gluon distribution / as a sum of contributions /„ from the production of n jets 

m = E Uv), (7) 

n=0 

then one gets a recursive relation 

f n (y)= [ y dy 1 K(y- yi )f n _ 1 (y 1 ), (8) 
Jo 

where K(y) is an y-dependent operator in the transverse momentum space 

K(y) = e yKuv K R . (9) 
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Eq. (8) allows to successively calculate the relative probabilities to produce n = 0, 1, 2, ... jets 
starting from the no-jet contribution determined by 

My) = e^-/(°)(0) + ( V dy^y-y^^-^. (10) 

Jo dyi 

In [4] the driving term was chosen to vanish at y = 0: 

f (0) (y,k) = A(l - e -yfe- k2 l Q » (11) 

(its normalization factor is irrelevant for our purpose). 

Distributions f n (y,k) themselves are not observable quantities. Physical probabilities 
are obtained by convoluting /„ with the gluon distribution in the projectile (the projectile 
impact factor). For the perturbative QCD to be applicable, a reasonable choice is to take the 
virtual photon as a projectile, as done in [4]. Having in mind that the BFKL picture may 
only be applied to low values of x, in our calculations we used a simplified expression for the 
virtual photon impact factor, independent of rapidity, which can be found in [5]. To have 
some qualitative idea of the situation in purely hadronic collisions, we have also made our 
calculations for a hadronic projectile with an unperturbative impact factor. For collisions of 
two identical hadrons it should be identical to the target impact factor which appears as an 
y-independent driving term f^°\k) in (2). 

In both cases the exclusive probabilities to observe n jets are given by 

_ f(d 2 k/k*)h(k)f n (y,k) 
^ n[y) f(d*k/k*)h(k)f(y,k)> {LZ) 

where h(k) is the impact factor of the projectile. Both impact factors, f^°\k) of the target 
and h(k) of the projectile, should vanish as k — > 0. This condition is satisfied by the virtual 
photon impact factor of [5]. As to the hadronic impact factor f(°\k), we have chosen it in 
close similarity with (11): 

f (0) {k)=k 2 e -kVQl (13) 

As with Eq. (11), the overall normalization is irrelevant. 

We are interested in the average values of < q > n in the observed jets, provided their 
number n is fixed. At this point one has to remember that the momentum k which serves 
as an argument of f(y, k) refers to the virtual gluon, and not to the emitted one, whose 
momentum q is hidden inside the kernel Kr. Therefore, to find an average of any quantity 
(j){q) depending on the emitted real jet momentum, one has to introduce the function (j)(q) 
into the integral (5), thus changing the kernel Kr to the kernel K av defined by 

(K av f)(k)=a s k 2 J ^^6(q 2 - ^^fih). (14) 

With n jets, one has to substitute one of the n operators Kr which generate the jets by K av , 
take a sum of all such substitutions, and divide by n. One has further to integrate over all 
momenta of the virtual gluon k multiplied by the projectile impact factor, and normalize the 
result to the total probability to have n jets. 

This recipe can be formalized in the following way. Introduce a generalized operator in 
the virtual gluon momentum space 

K 1 {y) = e yKuv [K R + K av \. (15) 

Let the function F(y, k) obey the equation 

F(y) = f (y) + T dyitfifo, - yi)F( Vl ). (16) 
Jo 
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One can split the function F into a sum of contributions F nm corresponding to the action of 
n operators K\, out of which m = 0, 1, ...n, are operators K a 



F(y) = E E F nm(y). (17) 

n=0 m=o 

Evidently F„o = f n - We are interested in the contribution F n \ = g n which contains a single 
operator K av . The average value of interest is determined by 

^, w l f(dk 2 /k*)h(k)g n (y,k) 

<m)>n~ nRdk2/kA)mUmky (18) 

In analogy with (8), one easily sets up a recursion relation for g n : 

g n = [ V dyiK(y - yi )g n -i(yi) + f dy^^^ K av f n ^( yi ), (19) 
Jo Jo 

with the initial condition go(y) = 0. Together with (8), this relation allows to calculate the 
function g n for n = 1,2, and then to use (18) to find the desired averages. 

The concrete choice of (f>(q) is restricted by the condition of convergence at large q: <f>(q) < 
q 2 , as q — > oo. For this reason we take <f>(q) = q, so that 

(iW) (k) = a s k 2 J ^r q e{q 2 - M 2 )/(&i)- (20) 

We defined our jets by taking \i = 2 GeV/c. As for the cutoffs, we used (3). To see the 
influence of the cutoffs we also repeated our calculations with Qq = 0.5 GeV/c and Qf = 1000 
GeV/c. The results slightly change in their absolute values (by no more than 1-6%) but both 
the n and y dependences remain the same. 

We have calculated the functions f n and g n from Eqs. (8) and (19) up to n = 5 and 
y = 15. Following [4] we have used the expansion in ./V Chebyshev polynomials to discretize 
the kernels in a simple way. The results we present have been obtained with TV = 80, although, 
as pointed out in [4], already iV = 20 gives a reasonable approximation. 

In Figs. 1-3 we present the averages < q > n for n = 1-5 and x = e~ y = 3.10~ 7 -0.1, for the 
7*-hadron collisions at Q 2 = 10, 100 and 1000 (GeV/c) 2 . In Fig 4 we show these averages 
for the collisions of two identical hadrons with the impact factor (13). As one observes, in all 
cases < q > n strongly grows with n at all rapidities. The growth is approximately linear 

<p T > n ^ a{y,Q 2 ) +b{y)n, (21) 

where both a and b depend on rapidity y, but, at fixed y, b is universal in the sense that it 
does not depend on Q 2 in DIS and has the same value for pure hadronic collisions. The slope 
b(y) falls with y: it is equal 1.25 GeV/c at y = 7.5, and 0.8 GeV/c at y = 15, so that one 
may expect that at ultrahigh energies < q > n will become independent of n. 

As an interesting byproduct of our study we find that the averages < q > n go down with 
rapidity for all n > 2. This is quite unexpected, since, as well-known, in the BFKL approach 
an overall average < q > rapidly grows with y ([6] and Eq. (25) below). It seems that this 
growth is totally explained by the growth of the average number of jets < n >. 

Passing to discussion, we first point out that it is an open question in which kinematical 
conditions and to what degree the BFKL pomeron may describe realistic hadronic processes. 
Emissions of high-pr jets in 7*-hadron collisions seem to be a suitable place to see the BFKL 
signatures. Our results show that in such emissions strong positive correlations are predicted 
between < px > and the number of jets, already for a single pomeron exchange. This in- 
dicates that in fact such correlations do not require multiple rescatterings nor pomeronic 
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interactions, but they are already present in the basic mechanism of jet production. Obvi- 
ously, this conclusion cannot be directly applied to particle production in the soft region, and 
so the empirically invoked absence of such correlations for particle production from the colour 
string does not really contradict our results. It can only be tested in its own framework by 
confronting colour string predictions with the experimental data. 

Our result (21) has been obtained for relatively small number of jets and at energies 
corresponding to y < 15. If extrapolated to all n and energies it would lead to a relation 
between the overall averages 

< p T >oc< n > . (22) 

However, it is well-known that this relation does not hold in the BFKL model at asymptotic 
energies, when < px > grows much faster than < n >. This more or less known fact can be 
demonstrated by a simple calculation. Indeed the inclusive cross-section I(y, y±, q) to produce 
a jet at rapidity y\ and with a transverse momentum q, in a collison with an overall rapidity 
y >> 1 is given by [7] 

/(y,i/i,9) = |^(v)^(l -*(*)), (23) 
where 



z = b\nq, b= -H -, a = 14a s C(3). (24) 

y ayi{y-yi) 

Here <3? is the error function and a(y) is the total cross-section. Since b ~ l/y 7 ?/ << 1, the 
term with <3? is actually important only at large Inq when it cuts the distribution in q at 
Inq ~ y/y. For this reason the scale of q is unimportant, so that one may safely fix it by 
setting fj, = 1. The average value of any positive power f3 of the transverse momentum is 
easily found to be 

z p ^ Jodyifd 2 qq' 3 I(y,y 1 ,q) 1 A 2 ,^ u x /-, /o\« / — roti\ 

< Pt >= -Vj — f , 2 j, r- = Tf e $(A), A = (1/2 /3^. 25 

Jo \ d Vi J d Q I {y,yi,Q) A 2 

Evidently < p T > grows exponentially with y for any f3 > 0. 

Thus our results cannot be valid for all n and energies and refer precisely to the values of 
n and energies for which the calculations were done. It is interesting to note that relations 
similar to (22), with < px > substituted by < p\ >, were earlier obtained from gluon 
saturation [8] , and in the percolation approach [9] . 

As mentioned, an unexpected result obtained in our calculation is that < p? > n at fixed 
n > 2 fall with energy. As seen in Figs. 1-4 this fall is not dramatic at energies at which 
we can expect the BFKL pomeron to be relevant (y > 10). Still, it is quite appreciable, 
especially in view of the naive belief that the average transverse momentum should grow 
with energy. This prediction can easily be tested experimentally as a possible signature of 
the BFKL pomeron. 

It is a pleasure to thank C. Pajares for interesting discussions and for pointing out to us 
some relevant references. We also thank J. Castro and G. Parente for their useful comments. 
This work was supported by the Secretarfa de Estado de Educacion y Universidades of Spain 
and by the RFFI grant 01-02-17137 (Russia). 
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Figure captions 

Figure 1. Average < pt > n for a fixed number n of jets produced in 7*-hadron collisions, as 
a function of x at Q 2 = 10 (GeV/c) 2 . Curves from bottom to top correspond to n = 1, 2, ...5. 
Figure 2. Average < px > n fo r a fixed number n of jets produced in 7*-hadron collisions, as a 
function of x at Q 2 = 100 (GeV/c) 2 . Curves from bottom to top correspond to n = 1, 2, ...5. 
Figure 3. Average < pt > n for a fixed number n of jets produced in 7*-hadron collisions, as a 
function of x at Q 2 = 1000 (GeV/c) 2 . Curves from bottom to top correspond to n = 1, 2, ...5. 
Figure 4. Average < px > n for a fixed number n of jets produced in hadronic collisions, as a 
function of y . Curves from bottom to top correspond to n = 1,2, ...5. 
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